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1 .  INTRODUCTION 

The  idea  of  bounds  on  the  optimal  value  of  a  geometric  program¬ 
ming  problem  has  been  present  since  the  inception  of  geometric  program¬ 
ming  (Duffin,  Peterson,  and  Zener  [5])  and  has  proved  to  be  fruitful  in 
applications  of  GP.  These  bounds  are  the  immediate  consequences  of  the 
duality  ..ueory  of  GP.  Recently,  this  idea  has  been  extended  to  paramet¬ 
ric  bounds  on  the  optimal  value  function  f*  .  Woolsey  (in  [2])  derives 
a  lower  bound  on  f*  using  the  known  fact  that  the  dual  objective  func¬ 
tion  at  a  fixed  dual-feasible  point  underestimates  f*  for  all  values 
of  coefficients  (parameters)  (see,  e.g.,  Dembo  [4]).  He  also  shows  how 
to  apply  this  result  in  a  practical  problem.  Fiacco  [7]  has  proposed  a 
general  approach  for  calculating  upper  and  lower  bounds  on  f*  (partic¬ 
ularly  simple  whenever  f*  is  convex  or  concave),  which  utilizes  sensi¬ 
tivity  information  as  well  as  Wolfe's  duality  theory.  This  paper  is 


T-464 


based  on  similar  Ideas,  except  that  GP  duality  theory  is  used  instead  of 
Wolfe's  duality,  and  the  special  structure  of  GP  primal  and  dual  problems 
is  exploited.  Several  classes  of  perturbed  GP  problems  are  shown  to 
possess  convex  or  concave  f*  or  at  least  "tight"  overestimating  and 
underestimating  problems  with  convex  f*  .  The  calculation  of  bounds  is 
illustrated  for  different  classes  of  perturbations  and  on  a  simple 
example  problem.  In  this  paper  we  are  mainly  concerned  with  posynomial 
GP  problems.  Possible  extensions  to  general  signomial  GP  problems  using, 
for  example,  the  idea  of  condensed  programs  (see,  e.g.,  [1])  remain  to  be 
developed.  Also,  the  topic  of  bounds  on  the  primal  and  dual  optimal  so¬ 
lution  points  is  not  discussed  here  (see,  e.g.,  [2,4]).  Bounds  on  f* 
based  on  a  general  idea  in  [7]  were  calculated  for  a  convex  equivalent 
of  a  GP  model  of  a  stream  water  pollution  abatement  system  by  Ghaemi  [11] 
and  by  Fiacco  and  Kyparisis  [9,10].  Bounds  on  f*  were  also  obtained  by 
Fiacco  and  Ghaemi  [8],  using  results  coinciding  with  some  results  of 
Section  4,  for  a  convex  equivalent  of  a  GP  model  of  a  power  system  energy 
model. 

2.  GENERAL  PRIMAL  AND  DUAL  BOUNDS  ON  THE  OPTIMAL 
VALUE  FUNCTION  OF  A  GP  PROBLEM 

A  posynomial  primal  geometric  programming  problem  is  a  nonlinear 

programming  problem  of  the  form 

min  gn(t,c) 
teEm  0 

subject  to  gk(t,c)  <1  ,  k=l,...,p  ,  PQ(c) 

tj  >  0  ,  j=l . m 


-  2  - 


where 
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m  a 


(t,c)  *  I  c.  II  t.1^  ,  k-0,l,...,p  , 

1  4—1  J 


l£Jk  ^ 


Dj^+1,  •  •  • ,  n^}  ,  k*0,l, . . .  ,p  , 


mQ  -  1,  «x  -  nQ+l,  m2  -  a^l,  mp  =  np_1+1»  np  =  n  » 

c  =  <clf...,cn)  . 

The  exponents  a.^  are  arbitrary  real  numbers  and  the  coefficients  c^ 

are  positive.  The  functions  are  posynomials. 

X1 

By  using  the  transformation  t^  =  e  J  we  derive  the  following 
equivalent  program 


rain  fn(x,c) 
xeEm  0 


subject  to  £,(x,c)  £1  ,  k=l,...,p 


P(c) 


where 


rm 

a . . x . 

f,  (x,c)  =  <  c.  e  ^  ^  ,  k=0,l,...,p 

K  i£j,  1 

k 

and  the  sets  are  the  same  as  in  Pq(c).  Program  P(c)  is  called  a 

transformed  primal  program  and  it  is  well  known  ([5])  that  it  is  a  con¬ 
vex  programming  problem  (for  any  fixed  c) . 

A  dual  geometric  programming  problem  corresponding  to  Pq(c)  (and 
P(c))  has  the  form 


-  3  - 
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c  ,  any  dual-feasible  point  <$  gives  us  the  lower  bound  on  f*(c)  and 
v*(c)  based  on  (1) ; 

f*(c)  _>  v*(c)  >  v(<5,c)  ,  V  c  .  (3) 

Woolsey  (in  Beightler  and  Phillips  [2])  utilizes  (3)  to  obtain  a  "tight” 
lower  bound  on  f*(c)  as  follows.  Suppose  that  is  the  optimal 

solution  of  D(Cq)  (it  is  necessarily  a  global  solution,  since  as  shown 
in  [5]  log  v(6,c)  is  concave  in  6  on  the  dual  feasible  set  for  any 
fixed  c  ).  Suppose  also  that  (A)  is  satisfied  so  that  (2)  holds.  Then 
from  (3)  and  (2)  we  obtain 

n  c.  ^Oi  n  c.  ^Oi 

f*(c)  >  v(5  ,c)  -  v(5  ,c  )  n  ~~  =  f*(c0)  n  ~  ,  V  c  .  (4) 

i*l  (C0ij  i=l  |_  Oi 

The  above  bounds  as  well  as  the  optimal  value  functions  f*(c)  and  v*(c) 
are  in  general  neither  convex  nor  concave.  In  the  next  sections  we  will 
identify  classes  of  problems  for  which  either  f*(c)  and  v*(c)  or  the 
bounds  on  f*(c)  can  be  shown  to  be  convex  or  concave.  This  is  impor¬ 
tant  since  it  simplifies  the  computation  of  the  bounds  on  f *(c)  (both 
upper  and  lower)  and  considerably  enhances  their  applicability.  These 
results  will  be  illustrated  using  the  following  example  problem  [5, 

p.  88]: 

—1  —1 

rain^  gQ(t,c)  -  ci  t^  *  t3  +  c^t.^  +  c3tit2t3 

subject  to  g1(t,c)  -  c^t"2  t“2  +  c5  t!f  t'1  £  1  EPQ(c) 

tt  >  0  ,  i=l ,2,3 
(cA  >  0  ,  1*1,2, ... ,5)  . 

-  5  - 
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The  transformed  problem  is  of  the  form 

-x.-hx~-x~  x.+x-  x  +x  -t-x 

min  fg(x,c)  =  c..  e  +  e  +  c.  e 

x£E3 


-2x1~2x2  JSx2-x3 

subject  to  f^(x,c)  =  e  +  e  £  1  . 


EP(c) 


The  dual  problem  corresponding  to  EP(c)  is 


max  v(S,c)  = 
6eE3 


' 

6, 

*  ' 

<S„ 

r  i 

6„ 

6. 

'  * 

C, 

1 

2 

3 

c , 

4 

Cc 

1 

2 

3 

4 

5 

6, 

5- 

6 . 

1J 

l  2J 

3J 

4J 

5j 

(«4  +  *5> 


<W 


subject  to  6^  +  62  +  <5^ 


"  61  +  62  +  63  "  2<S4 


=  1  , 
-  0  , 


EDO 


+  -  26^  +  ^6^  =  0  , 


"  61  +  62  +  63 


-  «5  =  0  , 


<Si  >  o  ,  i*l, 2, ...  ,5  . 


-1  !«  J-  2 

One  can  check  that  if  c^  c 2  c^  c^  >  2/27  ,  then  the  optimal  solution 


of  ED(c)  is  6*(c)  =  (1  -  2r(c),  r(c),  r(c),  2r(c)  -  h,  4r(c)  -  1)  , 

where  r(c)  = 


— ^  — V  —1  -1 

2  +  (4/27)c^  c2  2  c32  c^  c^  £  (*s, h)  and  the  optimal 


value  function  of  EP(c)  is  f*(c)  =  v*(c)  =  v(6*(c),c)  .  For  = 

(40,  20,  20,  1/3,  4/3)  ,  r(cg)  =  2/5  ,  the  optimal  dual  solution  is 
6*<c.)  =  (1/5,  2/5,  2/5,  3/10,  3/5)  ,  and  one  can  calculate  the  optimal 


primal  solution  x*(cn)  =  (0,  0,  log2) ,  and  f*(cQ)  =  v*(cn)  = 


'O' 


v(6*(Cq),Cq)  =  100.0.  The  lower  bound  (4)  is  thus  easily  applicable. 


-  6  - 
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3.  GP  PROBLEMS  WITH  CONVEX  OPTIMAL  VALUE  FUNCTION 

Suppose  that  coefficients  c^  are  functions  of  the  parameter 
vector  e  »  e  Er  of  the  form  c^(e)  =  l/h^(£)  ,  i-l,...,n  , 

where  h^(e)  are  positive  concave  functions  of  e  on  a  convex  set  Eq 
(if  h^(e)  =  h^p  ,  then  c^Ce)  *  1/h^Q  retna*-ns  unperturbed).  A  particu¬ 
larly  simple  class  is  obtained  when  h^(e)  =  >  0  ,  i=l,...,n  .  The 

primal  program  P(c(e))  can  be  written  as 

min  f.(x,e) 

xeEl“  P(£) 

subject  to  f^Cx.e)  ^  1  ,  k=l, — ,p 


where 


Vx-E>  'I  h 


ytn 

i  ,zj=i 


.  i  .X, 

jal  ij  j 


h— 0 ,l>.».»p 


iejk  -V 


If  we  denote  by  f*(e)  the  optimal  value  function  of  P(e),  then  the 
following  result  holds. 


Proposition  1.  ?*(e)  is  convex  on  E^  . 

Proof .  First  note  that  f^Cx.e)  can  be  written  as 

m 

f.(x,e)  =  l  exp  -log  h  (e)  +  l  a  ,x  ,  k=0,l,...,p  . 

K  ieJ.  j-1  J  J 

k  v  J 


Now,  introducing  new  variables  s^,»..,sn  ,  we  can  write  P(e)  in  an 
equivalent  form  as 


-  7  - 
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1 

i 

I 


in 


min  2  e  i 

(x,s)£Em*En  ieJQ 


r-  i 

subject  to  l  e  1 

ieJ, 


k-1, . . . ,p 


l  a . .x.  -  s.  <  log  h. (e)  ,  i=l,...,n  . 

j=l  2  2 

Programs  P(£)  and  P^(£)  have  the  same  optimal  value  function  f *(£)  . 

Since  h^(e)  is  concave  and  positive  for  £  e  E^  ,  log  h^Ce)  is  con-  . 

cave  for  e  e  E0  and  thus  P, (e)  is  a  jointly  convex  program  in  (x,s,£) 

J  i 

with  the  convex  optimal  function  f*(E)  for  £  £  [12].  Q.E.D. 


Proposition  1  enables  us  to  compute  the  upper  and  lower  bounds  on 
f*(e)  using  the  following  approach  proposed  by  Fiacco  [7],  Assume  that 
we  know  the  solutions  x^  of  P(eQ)  and  x.^^  of  P^),  Eq,  e.^  £  . 

Assume  also  that  the  conditions  of  the  sensitivity  theorem  (Fiacco  [6]) 
are  satisfied  at  e  =  £q  and  e  =  .  Then  the  gradients  of  f*(e) 

at  these  points  exist,  providing  us  with  the  lower  bounds  on  f*(e)  , 

E  e  [Eq.e^]  of  the  form  L^a)  =  f*^)  +  V£f*(Ei)  (e(ci)  -  £1>  ,  i=0,l  , 
where  £(a)  =  (1  -  ci)Eq  +  a£^  ,  a  E  [0,1]  .  The  upper  bound  on  f*(£(a)> 
is  given  by  U(a)  =  (1  -  a)f*(Eg)  +  af*(e^)  .  A  better  upper  bound  can 
be  obtained,  noting  that  f^(x,e)  are  jointly  convex  in  (x,e)  for 
e  £  Eq  .  This  implies  that  x(a)  =  (1  -  oi)Xq  +  ax^  is  feasible  for 
P(e(a)),  a  e  [0,1]  ,  so  that  f*(e(a))  fQ(x(a) ,E(a))  <  U(a)  ,  V  a  e  [0,1]  , 
the  last  inequality  following  from  convexity  of  fg(x(a) ,£(a))  .  It  is 
also  possible  to  derive  sharper  lower  bounds  on  ?*(e)  using  the  dual 
geometric  programming  problem  D(c(e)).  We  consider  this  problem  now. 


] 

I 

1 


-  8  - 
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The  dual  program  D(c(e))  has  the  form 

n  -6  P  X  (6) 

max  v(<5,e)  »  H  (h,(e)6  )  1  II  A  (6) 

6eEn  1«=1  k=l 

subject  to  X  (6)  -  l  5  ,  ic=-l,...,p  , 

K  ieJ, 

k 


l  1 

i£J0 


L  Vi  =  ° 


i=l 


5  >  0 


j=l, . . .  ,m  , 
i=l,...,n  . 


5(c) 


Denote  by  v*(e)  the  optimal  value  function  of  D(e).  Under  as¬ 
sumption  (A),  v*(e)  =  f*(e)  (provided  that  f*(e)  <  +°°)  and  this, 

together  with  Proposition  1,  implies  that  v*(e)  is  convex  on  E^  .  We 
will  prove  it  directly,  together  with  convexity  of  a  dual  lower  bound 
v(6,e)  ,  using  the  following  result. 


Proposition  2.  (i)  If  R  is  an  arbitrary  set  in  Em  ,  F(x,e)  is 

x 

concave  in  e  on  a  convex  set  S  C  E  for  any  fixed  x  e  R  ,  then 

F*(e)  =*  inf  F(x,e)  is  concave  on  S  . 

1  xeR 

(ii)  If  R  and  S  are  as  in  (i),  F(x,e)  is  convex 

in  E  on  S  for  any  fixed  x  £  R  ,  then  F*(e)  =  sup  F(x,e)  is  convex 

xeR 


on  S  . 


Proof .  (i)  Let  e^^  e  s  »  X  e  [0,1]  .  Then 


F*(X£l  +  (1  -  X)e2) 


inf  F(x,  Xe..  +  (1  -  X)e~) 
xeR 


>  inf  XF(x,e1)  +  (1  -  A)F(x,e„) 

^  xeR  1 

X  inf  F(x,e-)  +  (1  -  X)  inf  F(x,e2) 
xeR  xeR 


«  X  FJ(ex)  +  (1  -  X)F*(e2)  . 
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(ii)  F*(e)  =  sup  F(x,£)  =  -inf (-F(x,£))  .  Since  -F(x,e) 
xeR  xER 

satisfies  assumptions  of  (i),  -F*(e)  is  concave  on  S  by  (i),  so 

F*(e)  is  convex  on  S  .  Q.E.D. 

Remark.  Proposition  2  appears  to  be  a  variation  of  a  well  known  result 

of  convex  analysis:  if  all  4>  e  ♦  are  convex,  <J>  is  an  arbitrary  set, 

then  4>(x)  =  sup  <|>(x)  is  also  convex  (on  a  convex  set). 

<J>e$ 


Proposition  3.  (i)  v(5,e)  is  convex  in  e  on  for  any  fixed 

dual-feasible  6  . 


(ii)  v*(e)  is  convex  on  E^ 
Proof .  (i)  v(S,e)  can  be  written  as 


v(<5,e)  =  exp 


[-  jL  Sl 


n  a,«) 

k=l  * 


Xk(6) 


(We  define  x  log  x  *  0  for  x  *  0  .)  Since  h^e)  is  concave  and 

positive  on  Eq  ,  -6^  log(h^(e)6Q)  is  convex  on  E^  for  any  dual- 

feasible  6  and  thus  v(6,e)  is  convex  in  e  on  E^  for  any  such  6  . 

(ii)  Denote  the  dual  feasible  set  by  R^  .  Then  v*(e)  =  sup  v(6,ei) 

6£RD 

J  _  _ _ ___ _  _ _ .  n  1 _  /  »  \  _ t  r» _ _  t  .  r  <-»  /  <  »  \ 


is  convex  on  Eq  by  (i)  and  Proposition  2(ii) 


Q.E.D. 


Suppose  now  that  condition  (A)  is  satisfied  so  that  f*(e)  =  v*(e)  . 
Let  6q  and  6^  be  the  solutions  of  D(Eq)  and  D(e^),  respectively.  If 
the  assumptions  of  the  sensitivity  theorem  hold  for  either  problem  P(e) 
or  D(e)  at  £  =  £q  and  £  =  £^  ,  then  the  gradients  of  both  f*(e)  and 
v*(e)  exist  at  £  *  Eq.e^  and  we  have 

V£***eiJ  =  Ve^*(£i)  =  Ve^(6i,ci)  »  1=0,1  '  (?) 
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Since  also  f*^)  -  v*^)  -  ,  i»0,l  ,  convexity  of  v^.r.)  , 

i=0,l  implies  that  for  i»0,l  ,  a  £  [0,1]  , 

Lt(a)  -  v(6i,ei)  +  Vev(5i,e1)(e(a)  -  c±)  <  v(6i,e(a))  ,  (8) 

proving  that  v(6^,e(a))  ,  i=0,l  ,  are  uniformly  better  lower  bounds  on 
f*(e(a))  than  L^Ccx)  .  In  fact,  to  derive  the  above  bounds  we  only  need 
the  directional  derivatives  of  f*(£)  =  v*(£)  in  directions 
and  £q  -  £  .  If  the  feasible  set  of  D(£)  is  compact  and  unique  global 

solutions  Sq  of  D(£q)  and  <5^  of  D(e^)  exist,  then  from  Danskin's 
theorem  [3]  we  obtain  the  directional  derivatives  of  f*(£)  at  e  = 
in  the  direction  z  as 


Dzf*(£1)  -  Dzv*(£i)  =  Vev(6i,£i)z  , 


i=0, 1  . 


(9) 


Consider  now  a  more  general  class  of  perturbations  of  the  form 


c±(e) 


1  6i£ 

n  h.J*(e) 
£==1  * 


_1  R 

1  ~^i£ 

-  n  h  (e) 

i-i 


where  h^(e)  are  positive  concave  functions  of  e  on  a  convex  set 


E0CE 


,  £^  are  positive  integers,  and  8^  >_  0  for  £=1,...,£^  , 

,n  .  (The  previous  case  is  obtained  by  setting  £^  =  1  ,  8^  =  1  , 


i=l,...,n  .)  All  results  obtained  in  this  section  extend  immediately  to 

*i 

this  class  (note  that  log  c^(e)  *  log  h^(£)  convex  on 

Eq  ).  This  allows  us  in  particular  to  obtain  bounds  on  convex  f*(e) 

r  -^i£ 

for  perturbations  of  the  type  c^e)  =  ,  8^  0  ,  £=l,.,.,r  , 

i*l, . . . ,n  , 

Oembo  [4]  considers  a  slightly  more  general  problem  than  P(c)  of 
the  form 


min  fn(x,c) 
xeE”1  0 

subject  to  f^(x,c)  <  r^  ,  k*l,...,p 


P(c,r) 
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where  ,  k=l,...,p  are  positive  numbers  and  functions  , 
k=0,l,...,p  are  defined  in  P(c).  If  we  define  c^(e,r)  =  c^O/r^  , 
i  e  ,  k=l,...,p  ,  where  c^e)  has  the  general  form  given  above, 
then  it  is  easy  to  show  that  all  results  of  this  section  extend  to  this 
problem,  too.  This  is  also  clear  from  the  formulation  of  P(c,r)  and  its 
corresponding  dual  D(c,r)  in  terms  of  (c,r)  , 

n  (c.)6i  P  [A.(6)}V6) 

max  v(6,c,r)  =  II  -g—  H  -  D(c,r) 

SeRp  1=1  (  ij  k=l  rk 

where  is  the  feasible  set  of  D(c)  and  are  defined  in  D(c)  . 

As  an  illustration  of  the  preceding  results,  consider  our  example 
problem,  EP(c),  with  the  following  perturbations:  c(a)  =  (c^/ (1+a) , 
c02>  c03»  co4*  CQc/d+a))  >  a  G  [0,1]  •  The  optimal  value  function 
f*(a)  of  EP(c(a))  is  convex  on  [0,1]  by  Proposition  1.  Since  c(0)  = 
cQ  ,  f *(0)  =  v*(0)  *  100.0  ,  6Q  =  6*(c(0))  =  (1/5,  2/5,  2/5,  3/10,  3/5) 
and  xQ  «  x*(c(0))  =  (0,  0,  log2)  ,  where  v*(a)  denotes  the  optimal 
value  function  of  ED(c(a)).  Since  r(c(a))  =  2/(5+a)  ,  for  a  =  1  , 
c1  =  c(l)  =  (20,  20,  20,  1/3,  2/3)  ,  r(c(l))  =  1/3  ,  and  =  6*(c(l)) 

(1/3,  1/3,  1/3,  1/6,  1/3)  .  Calculations  give  us  x^  =  x*(c(l))  = 

(0,  0,  0)  and  f*(l)  *  v*(l)  =  v(A^,c^)  =  60.0  .  Define  v(6,a)  = 
v(fi,c(a))  .  From  inequality  (4)  we  obtain  convex  dual  lower  bounds 

,  z.,„, f  1  )S01+S05  100.0 

'*<“>  i  ’ E  <0)lm;j  ■  -  -  o).8 

and 


f*(a)  >  v(61,a)  =  f*(l)  ~ — 


611+615  22/3  *  60.0 


(1  +  a) 


2/3  ’ 


a  c  [0,1]  . 


The  linear  lower  bounds  on  f*(a)  are  computed  using  the  formulas  for 
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v(61,a)  ,  i-0,1  ,  as 

LQ(a)  =  f*(0)  +  \7av(60,0)a  -  100.0  -  80.0a 

and 

Lx(a)  -  f*(l)  +  Vav(61,l)(a-1)  -  60.0  -  20.0(a-l)  =  80.0  -  20.0a  . 

The  linear  upper  bound  is  computed  as 
U(a)  -  (l-a)f*(0)  +  af*(l)  -  (l-a)lOO.O  +  a60.0  =  100.0  -  40. 0a  . 

A  sharper  convex  upper  bound  is  given  by 

fQ(a)  -  fQ(x(a) ,c(a))  , 


where  x(ct)  =  (l-a)xg  +  ax^  ,  a  e  [0,1]  and  calculations  show  that 
fg(a)  *  40.0(1/2)*  a/(l+a)  +  40.0*2*  a  .  All  the  above  bounds  are  de¬ 
picted  in  Figure  1. 
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Consider  now  another  class  of  perturbations  which  yields  a  convex 
optimal  value  function  f*(c)  •  Assume  that  only  the  coefficients  of  the 


objective  function  f^x.c)  ,  c  =  (c^,...,cnQ)  ,  are  perturbed  and  define 
ci(e)  =  Cq-jY^e)  ,  3i  >  1  »  i=l,...,nQ  ,  where  y(e)  is  a  positive  con¬ 
vex  function  of  e  on  a  convex  set  C  Er  .  Note  that  all  the  coeffi¬ 
cients  in  fg(x,c)  are  perturbed  now. 


Proposition  4. 

(i)  For  any  fixed  primal-feasible  x  ,  fg(x,c(e))  is  convex  in  e  e 

(ii)  For  any  fixed  dual-feasible  6  and  fixed  c  =  (c  . c  )  , 
v(e)  =  v(S,c(e),c)  is  convex  in  e  e  EQ  . 

(iii)  The  optimal  value  function  of  D(c(e))  ,  v*(e)  ,  is  convex  in 
e  c  Eg  and  if  condition  (A)  also  holds,  the  optimal  value  function  of 
P(c(e))  ,  ?*(e)  ,  is  convex  in  e  e  Eg  . 

Proof. 

(i)  Follows  immediately  from  the  convexity  of  y(e)^  for  any  g  >  1 
and  the  form  of  fg(x,c)  . 

(ii)  Since 


'Ke) 


n0 

n 

i=l 


Si' 

5i 

C01Y<£) 

n 

n 

i3n0+1 

C  . 

1 

6i 

> 

l«J 

P  Ak(6) 

n  A.  (6) 

k=l  K 


we  can  write 


00  Bi5i  k~l  Bi6i 

v(e)  =  A  It  y(e)  1  1  =  Ay(e)  11  11 

i=l 


(10) 


where 


“0 

n 

i=l 


O 

O 

i  n 

n 

' 

ci 

ITT 

T~ 

{  iJ 

i=no+1 

ij 

P 

n 

k«l 


\(«) 

Ar(6) 
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n  nQ  nQ 

Since  6  is  dual-feasible,  J  ^  5^^  *=  1  ,  so  that  £i=1  6^  >  £.=1  61  =  1 

6 

and  the  result  follows  from  the  convexity  of  Y(e)  ,  3  >  1  . 

(iii)  The  first  assertion  of  this  part  follows  immediately  from  part 
(ii).  Proposition  2(ii),  and  the  fact  that  the  dual  feasible  set  R^  is 
fixed.  The  last  assertion  is  a  consequence  of  the  equality  f*(e)  *  v*(e) 
under  condition  (A).  Q.E.D. 

This  result  enables  us  to  calculate  bounds  on  f*(e)  in  a 
similar  way  as  in  the  first  part  of  this  section.  Using  Danskin's 
theorem,  one  can  find  the  directional  derivatives  of  f*(e)  using  the 
gradients  of  v(e)  and  obtain  linear  lower  bounds  on  f*(e)  .  The 
linear  upper  bounds  U(e)  can  also  be  obtained  as  before.  Sharper 
convex  bounds  v(e)  are  available,  too.  However,  since  fg(x,c)  is  in 
general  not  jointly  convex  in  (x,c)  ,  the  upper  bound  fQ(x(a) ,c(e(a))  , 
x(a)  =  (l-a)xQ  +  ax^  ,  e(a)  =  (l-a)eQ  +  ae^  ,  a  e  [0,1]  will  not  neces¬ 
sarily  be  convex  and  better  than  the  linear  upper  bound.  The  convex 
upper  bounds  fgCx^.cCefo)))  ,  i“0,l  ,  will  be  better  than  U(e(a))  only 
for  a  close  to  0  (i»0),  or  1  (i*l)  in  general. 

This  approach  can  be  extended  to  include  perturbations  in  the  coef¬ 
ficients  c  of  the  constraint  functions,  also  of  the  form  c . (e)  = 

Bi 

CQiY(e)  ,  but  with  8^  ^  0  ,  i  *  n^+l . n  (note  that  8^  1  for 

i*l, . . . ,ng) .  It  can  be  easily  shown  that  Proposition  4  remains  true  in 
this  case.  However,  since  the  primal  feasible  set  now  depends  on 

e  ,  at  least  one  of  the  solutions  x^  of  P(c(Eq))  or  x^  of  P(c(e^)) 
will  no  longer  be  feasible  in  general  for  all  values  of  e  e  [eQie^]  , 
reducing  the  availability  of  upper  bounds  f^Cx^cCe))  ,  i*=0,l  .  Also, 
in  general  x(a)  will  not  be  feasible,  so  that  f^xCa ),c(e(a)))  cannot 
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be  used  as  an  upper  bound  on  f*(£)  .  The  linear  bounds  on  f*(t)  and 
the  dual  lower  bounds  will  nevertheless  remain  valid. 

Even  more  generally,  Proposition  4  will  be  valid  whenever 
£$^<S*(c(e))  >  1  for  the  considered  perturbations.  We  utilize  this 
fact  in  the  following  example  where  it  can  be  shown  that  =  1  , 

a  e  [0,11  . 

We  now  consider  the  example  EP(c),  with  perturbations  c(a)  = 

(cQ1(l  -  a/2),  cQ2,  cQ3(l  -  a/2)2,  c^(l  -  a/2),  cQ5)  ,  a  e  [0,1]  .  The 
optimal  value  function  f*(e)  of  EP(c(a))  is  convex  on  [0,1]  by  Proposition 
4(iii)  (we  set  y(a)  -  1  -  a/2,  ^  =  1,  B2  -  0,  B3  -  2,  B4  =  1,  B5  =  0)  . 

Since  c(0)  «  cQ,  f*(0)  =  v*(0)  =  100.0,  6Q  =  <5*(e(0))  =  (1/5,  2/5,  2/5,  3/10, 
3/5)  ,  and  Xq  =  x*(c(0))  *  (0,  0,  log2)  .  Since  r(c(a)>  =  (2-a)/(5-2a)  , 
for  a  =  1  ,  c^  =  c(l)  *  (20,  20,  5,  1/6,  4/3)  ,  r(c(l))  =  1/3  ,  and 
6^  =  6*(c(l))  *  (1/3,  1/3,  1/3,  1/6,  1/3)  .  We  calculate  =  x*(c(l))  = 
(log(1.05),  log4,  log (.67))  and  f*(l)  -  v*(l)  =  v(6  Cj_)  =  42.43  . 

Inequality  (4)  yields  the  following  convex  dual  lower  bounds  for 
a  c  [0,1]  , 


f*(a)  > 


v(6Q,a) 


f*(0) 


M 


601+2<503+604 


100.0 


f) 


1.3 


and 


un  uiruu  7/a 

f*(a)  >  v(6va)  =  f*(l) (2  -  a)  11  =  42.43(2  -  a)  '  . 

The  linear  lower  bounds  on  f*(a)  are  calculated  with  the  help  of 
v(6i,a)  ,  i=»0,l  ,  as 

LQ(a)  =  f*(0)  +■  Vav(60,0)a  «  100.0  -  65.0a 
and 
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Lx(a)  -  f*(l)  +  V  v(6ltl)(a  -  1)  *  42.43  -  22.05(a-l) 

±  64.48  -  22.05a  . 


The  linear  upper  bound  Is 

U(a)  -  (l-a)f*(0)  +  af*(l)  -  (l-a)lOO.O  +  a42.43 

A  100.0  -  57.57a  . 

Since  c^a)  are  nonlncreaalng  in  a  ,  the  primal  feasible  set  is 
increasing  in  a  .  Thus  the  optimal  solution  point  xQ  of  P(c(0))  is 
primal-feasible  for  all  a  e  [0,1]  .  Since  the  constraint  f^x.c)  is 
binding  at  the  optimal  solution  point  of  P(c(l))  for  0  =  0^,  x^ 

is  not  primal- feasible  for  any  value  of  a  e  [0,1)  .  Therefore,  we  can 
only  use  the  convex  upper  bound 

f0(x0,c(a))  -  20.o|l  -  |]  +  40.0  +  40.o[l  -  §]2 

-  10. 0a2  -  50. 0a  +  100.0  . 

We  depict  the  above  bounds  in  Figure  2  (f*(a)  is  not  depicted,  since  the 
bounds  U(a)  and  v^.a)  ,  i»0,l  are  very  tight). 


4.  GP  PROBLEMS  WITH  CONCAVE  OPTIMAL  VALUE  FUNCTION 

Consider  again  the  problem  P(c)  and  assume  that  we  perturb  the 

coefficients  c^  in  the  objective  function  fg(x,c)  only.  Denote  the 

vector  (c. , • • . , c  )  by  c  ,  the  optimal  value  function  of  P(c),  by 
i.  no 

f*(c)  and  the  fixed  feasible  set  of  P(c)  by  Rq  . 

Proposition  5.  f*(c)  is  concave. 

Proof.  Follows  immediately  from  Proposition  2(i)  since  fg(x,c)  is 
linear  in  c  for  any  fixed  x  .  Q.E.D. 
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This  result  enables  us  to  derive  the  upper  and  lower  bounds  on 
f*(c)  ,  given  the  solutions  Xq  of  P(Cq)  and  of  P(c^)  and  the 

gradients  of  f*(c)  at  c  =  Cq  and  c  =  ,  in  the  same  manner  as  it 

was  done  in  Section  3.  However,  the  upper  bound  fg(x(a) ,c(a))  , 
x(a)  =  (1  -  oOxq  +  otx^  ,  c(a)  =  (1  -  <x)cq  +  ac^  ,  a  e  [0,1]  ,  will  be 
neither  concave  nor  convex  in  general  and  will  not  necessarily 


Legend 


Figure  2. — Bounds  on  the  convex  f*(a) 
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underestimate  U^(a)  .  As  mentioned  before,  directional  derivatives  of 
f*(c)  at  c  ■  co,cl  are  8uf f icient  for  our  purposes.  Thus,  if  Rq  is 
compact  and  x^  and  are  the  unique  global  solutions  of  P(Cq)  and 

P(c^)  ,  respectively,  then  by  Danskin's  theorem  [3]  we  obtain  that 

V*^)  -  V-fQ(xi,ci)z  =  f0(x±,z)  ,  i=0,l  (11) 

which  means  that  the  upper  bounds  on  f*(c(a))  are  of  the  form  IL(a)  = 

f()(Xi,c(&))  ,  i=0,l  ,  a  £  [0,1]  .  If  the  condition  (A)  holds,  then  the 

optimal  value  function  of  D (c)  ,  v*(c)  ,  is  equal  to  f*(c)  and  is  thus 

concave  by  Proposition  5.  Now  we  prove  concavity  of  the  lower  bound 

v(6,c,c)  in  c  for  any  6  and  c  =  (c  c  )  fixed. 

n0+1  n 

Proposition  6.  For  any  fixed  dual-feasible  6  and  fixed  c  ,  v(<S,c,c) 

is  concave  in  c  . 

n0  6i 

Proof.  Denote  vq(c)  *  v(6,c,c)  ■  A(6,c)  11^^  c i  ,  where 

n  6  n  -5  P  X  (6) 

A(6,g)  *  n  c  1  n  6  1  n  \  (6)  K  >  o 

i«n0+l  1  1=1  1  k=l 


is  fixed.  Calculations  give  us  V-v_(c)  =  [6,/c,  ,  . 6  /c  ]v_.(c) 

c  u  i  I  n0  nQ  u 


V-Vq(c)  =  Vq(c)  *  v  (c) 


VC1 


6  /c 

.  n0  n0. 


[6  /c  ,  ....  5  /c  ]  - 
11  n0  n0 


'Vcl 


0  6  /t 


n0  n0 


(12) 


T  2  -  un  - 

We  want  to  show  that  y  V-Vg(c)y  <  0  ,  V  y  s  E  .  Since  Vq(c)  >  0  ,  it 

it  enough  to  show  that 
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6  /c 

r  l 


6  /< 


(6  /c  ,  ....  6  /c  ]y  -  y 
11  no  n0 


L  no 

noJ 

n0  6  y . 
c  i7i 

2  n.  A 

0  o  .y 
V  1 

Ci, 

4  2 

i=l  cT 
i 

Yci 


6  /c2 

n0  n0J 


<  0  ,  V  y  e  E 


(13) 


But  this  follows  immediately  from  the  Cauchy-Schwarz  inequality  [5] 


n 

r 

2 

f 

r  2 

[  "  ,2 

/  n.  b. 

< 

/  a 

y  b. 

k=l  k  k 

L  ab 

[k=i  kJ 

[k-l  J 

V  a,b  e  E 


(14) 


if  we  set  n  =  ,  a^  =  6^  ,  =  6^  (y^/c^)  >  k=l,...,n  and  recall 

that  Ik*i  <$k  *  IkeJ^  \  “  1  •  Q.E.D. 


If  6^  uniquely  solves  D(Cq),  6^  uniquely  solves  and  if 

the  dual  feasible  set  is  compact,  then  under  the  assumption  (A), 

f*(c)  =  v*(c)  ,  and  thus  by  Danskin's  theorem  for  any  z  , 

Dzf*^ci^  =  Dzv*^ci^  =  VcV^i’Ci’^Z  ’  i=0,l  (13) 

This  shows  that  the  lower  bounds  v(6^ ,c (a) ,c)  ,  i=0,l  ,  are  at  least 
initially  better  than  the  linear  lower  bound  L(a)  =  (1  -  a)f*(c^)  + 
af*(Cj)  but  in  general  they  do  not  have  to  be  uniformly  better  for  all 
a  G  [0,1]  . 

There  are  two  ways  to  extend  the  results  of  this  section.  One  way 
is  to  assume  that  c^(g)  =  h^(e)  ,  i=l,...,ng  ,  where  h^(e)  are  concave 
positive  functions  of  the  parameter  vector  e  =  (e^,...,Gr)  .  The  case 
treated  before  is  obtained  when  we  set  h^(c)  =  ,  i=l,...,nQ  .  One 

can  show  that  Propositions  5  and  6  continue  to  be  true  for  this  more 
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general  case.  Upper  and  lower  linear  bounds  can  be  obtained  as  before, 
given  the  solutions  and  directional  derivatives  of  the  optimal  value 
function  ?*(e)  at  e  -  eo»ei  *  objective  function 

n.  cm 

ffl(x,e)  =  l  h  (e)  e  J=  1J  J 
i-1 

gives  us  now  a  sharper  upper  bound  than  the  linear  one,  since  from  con¬ 
cavity  of  f^Cx.e)  in  £  for  any  fixed  x  ,  fixedness  of  R  ,  and 
equality  Dzl*(e^)  ■  V^f^x^.E^z  »  i-0,1  ,  it  follows  that 

f*(e(a))  <  f0(xlte(a))  <  Ut(a)  ,  i-0,1  (16) 

where  e(a)  *  (1  -  a)£g  +  ae^  ,  a  e  10,1]  and 

Ui(a)  -  f*(ei)  +  7ef0(xi,e1)(e(a)  -  ,  i-0,1  . 

The  dual  lower  bound,  concave  by  the  extension  of  Proposition  6,  applies 
as  before.  Another  generalization  is  obtained  by  noting  that  Proposition 
5  is  valid  even  if  some  coefficients  c^  are  negative  and  directions  of 
inequalities  defining  the  feasible  set  changed.  In  other  words,  the  op¬ 
timal  value  function  f*(c)  is  concave  also  in  case  of  the  general  sig- 
nomial  geometric  program.  In  this  case  we  can  still  utilize  linear  bounds 
on  f*(c)  based  on  sensitivity  information  [7],  although  the  dual  bounds 
may  no  longer  be  valid.  In  order  to  illustrate  the  above  bounds,  we  use 
the  same  example  EP(c)  as  before. 

Define  the  perturbations  as  follows:  c(a)  -  (Cq^(1  -  a/2), 

Cq2(1  -  3a/4),  Cqj,  Cq^,  Cq^)  ,  a  £  [0,1]  .  Since  c(a)  is  linear  in  a 
this  is  the  basic  case  and  the  results  of  this  section  apply,  so  that  the 
optimal  value  function  f*(a)  of  EP(c(a))  is  concave  on  [0,1]  by 
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Proposition  4.  Since  c(0)  =  Cq  ,  f*(0)  -  100.0  ,  -  (1/5,  2/5,  2/5, 

3/10,  3/5)  ,  and  xQ  =  (0,  0,  log2)  ,  as  determined  before.  Since 
r(c(a))  *  (2  +  (1  -  a/2)//4-3a)_1  ,  for  a  =  1  ,  cJ  =  c(l)  *  (20,  5,  20, 
1/3,  4/3)  ,  r (c (1 ) )  *  2/5  ,  and  ^  “  6* (c (1) )  =  6Q  .  Calculations  give 
Xi  =  x*(c(D)  »  (log4 ,  -log4,  0)  and  f*(l)  =  v*(D  =  v(61,c]L)  =  50.0  . 
Define  v(6,a)  =  v(6,c(a))  as  before.  Inequality  (4)  provides  concave 
dual  lower  bounds  for  a  £  [0,1]  , 

f*(a)  >  v(6Q,a)  =  f*(0)[l  -  |  01(l  -  ^  02  =  100. o[l  -  |  1  -  ^ 


and 


f*(a)  >  v(6lta)  =  f*(l) 


fi  -  f]H  -  ff02 


=  100.0  1  - 


(note  that  v(6Q,a)  =  v^.a)  since  6Q  =  6^.  The  linear  lower  bound 
is  given  by 

L(a)  =  (1-a) f*(0)  +  af*(l)  =  (l-a)lOO.O  +  a50.0  =  50.0  -  50.0a  . 

The  linear  upper  bounds  are  given  by 

UQ (a)  =  fQ(x0,H(a))  =  100.0  -  40.0a 

and 

U1(a)  =  fQ(x^,c(a))  =  120.0  -  70.0a  . 

We  can  also  calculate  the  upper  bound, 

fQ(x(a),c(a))  =  20.0 |l  -  | ]  +  20.o[l  -  —■  21+a  +  20.0-21_a  . 

We  depict  the  calculated  bounds  in  Figure  3.  (The  upper  bound 
fp(x(a) ,c(a))  is  not  depicted  since  it  is  almost  equal  to  f*(a)  .) 
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5.  GP  PROBLEMS  HAVING  JOINTLY  CONVEX  OVERESTIMATING 
AND  UNDERESTIMATING  PROBLEMS 

Consider  now  frequently  used  linear  perturbations  of  the  coeffi¬ 
cients  c^  in  both  the  objective  function  as  well  as  in  the  constraint 
functions  of  the  program  P(c).  This  "natural"  class  of  perturbations  is 
not  covered  by  any  case  considered  so  far.  The  optimal  value  function 
f*(c)  is  in  general  neither  convex  nor  concave  in  this  case.  However, 
it  is  possible  to  define  "tight"  overestimating  and  underestimating  prob¬ 
lems  with  convex  optimal  value  functions  to  which  in  turn  bounding  tech¬ 
niques  of  Section  3  can  be  successfully  applied.  Specifically,  consider 

perturbations  of  the  form  c.(a)  =  c„.(l  -  d.a)  ,  where  d.  >  0  , 

iOii  i  = 

i=l,...,n  and  a  £  [0,a  ]  for  some  0  <  a,.  <  min  (1/d.)  (we  define 

l£i£n  1 

1/d^  =  +oo  if  d^  =  0  ) .  Obviously,  if  d^  =  0  for  some  j  ,  then  the 
coefficient  c^  will  not  be  perturbed.  Using  linear  bounds  on  the  con¬ 
vex  function  1/(1  -  d^a)  on  the  interval  [0,ctg]  ,  one  can  easily  de¬ 
rive  the  following  inequalities: 


max  1  +  d^a. 


1  +  d.(a-2a0) 

U-Vq)2 


1  +  d  (a-a  ) 

<  - ± -  <  - 1 - 

==  1  -  d.a  =  1  -  d,a„ 

i  i  0 


and  (substitute  +°°  for  (l-d^a^)  / (1  +  d^(a-2aQ>)  when 
2 

a  £  a  max  (d . )  ) 
l<i<n  1 


C0i(1~di°t0)  ,  ,  ,  .  c0i  c0i(1  dia0)  . 

1  +  d1(a-aQ)  =  ci^a  ~  c0i  1-dia^  =  min  i  +  d^’  1  +  di(a-2aQ)  (18) 


Define  the  following  geometric  programming  problems: 


min  fQ1(x,a)  s.t.  fj^x.a)  <  1  , 
x£Em 


k=l , . . . , p 


PU1(ot) 
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From  (18)  it  follows  that  for  a  e  [0,a^]  ,  x  e  Em  , 

f£(x,a)  <  fk(x,c(a))  <  minjf^Cx.a) ,  f^2(x,a)|  ,  k=0,l,...,p  (22) 

This  inequality  asserts  that  both  PU^(a)  and  PU2(a)  are  the  overestimat¬ 
ing  problems  for  P(c(a))  and  that  PL(a)  is  the  underestimating  problem 
for  P(c(a)),  where  a  e  [0,ag]  . 

Denote  by  fj^(a)  ,  f*^(a)  ,  and  f*^(a)  the  optimal  value  func¬ 
tions  of  the  programs  PU^(a),  PU2(a),  and  PL(a),  respectively.  (Note 
that  f£U(a)  ■  for  a  e  [0,5q]  . )  From  inequality  (22)  we  obtain  for 
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a  e  (0,aQ] 

f*L(a)  <  f*(c(a))  <  minjf*U(ct),  f*U(ct)J  (23) 

Since  f^(x,0)  =  f^(x,c(0))  =  f^(x,0)  and  f^(x,a^)  =  f^(x,c(aQ))  = 
f^2(x,aQ)  for  k=0,l,...,p  and  x  e  Em  ,  we  also  have  that 

f*L(0)  =  f*(c(0))  =  f*U(0)  (24) 

and 

f*L(aQ)  =  f*(c(aQ))  =  f*VQ)  (25) 

which  proves  that  the  bounds  on  f*(c(a))  are  "tight"  at  a  =  0  and 
a  =  a0  . 

The  appealing  feature  of  problems  PU^(a),  Pl^a),  and  PL(a)  is 
that  they  are  jointly  convex  in  (x,a)  (Pl^a)  is  jointly  convex  in  the 
extended  sense)  and  therefore  possess  convex  optimal  value  functions 
f*U(a)  ,  f*U(a)  ,  and  f*L(a)  ,  respectively.  This  is  immediately  seen 
from  the  fact  that  the  coefficients  in  all  the  above  problems'  functions 
are  of  the  form  l/h^a)  ,  where  h^(a)  are  linear  positive  functions  of 
a  ,  and  Proposition  1  of  Section  3.  Since  f*^(a)  underestimates 
f*(c(a))  ,  convex  dual  lower  bounds  as  well  as  linear  lower  bounds  on 
f*L(a)  can  be  calculated  using  results  of  Section  3  and  then  used  as 
lower  bounds  on  f*(c(a))  .  Similarly,  convex  or  linear  upper  bounds  on 
f*^(cx)  and  f^U(a)  can  be  computed  also,  using  results  of  that  section 
(in  the  case  of  f*^(a)  for  a  such  that  f*^(ot)  <  +°°  ) ,  providing  us 
with  upper  bounds  on  f*(c(a))  by  virtue  of  inequality  (23). 

The  approach  described  here  can  be  extended  to  the  case  where  the 
c^'s  are  concave  and,  e.g.,  decreasing  functions  of  a  parameter  a  .  We 
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will  not  discuss  this  idea  in  detail  here,  since  the  derivation  of  bounds 
in  this  case  is  similar  to  one  described  above. 

Consider  once  more  our  example  problem  EP(c),  with  perturbations 

c(a)  *  (cQ1(l  “  (5/8)a),  Cq2»  cq3»  c04^  ~  a/4),  c05^1  "  ct'/2^  *  Since 

d  «  5/8  ,  d  •  0  ,  d  =  0  ,  d  **  1/4  ,  d  =  1/2  ,  min  (1/d.)  =  8/5  , 

J  4  J  l<i<5  1 

and  we  choose  *»  1  .  Since 


r(c(a)) 


2  +  - 
2 


>-f  ‘-f 


-1 


for  a  =  1 


cx  =  c(l)  -  (15,  20,  20,  1/4,  2/3)  ,  r(c(l))  =  1/3  ,  and  =  5*(c(l))  = 

(1/3,  1/3,  1/3,  1/6,  1/3)  .  From  calculations  we  obtain  =  x*(c(l))  = 

(log(.86),  0,  0)  and  f*(l)  ■  v*(l)  *  v(6^,c1)  =  51.96  .  The  overesti¬ 
mating  problems  are  of  the  form 

it  c0l  **xl”^JC2~X3  Xl+X3  xi+x2+x3 

min  fj  (x.a) - e  1  2  3  +  c._  e  1  3  +  c..  e  1  2  3 

xeE3  l  +  |  a  UJ 


U  c04  -2xi-2x2  c05  }5x2~x3 

subject  to  f”,  (x,o)  -  e  1  2  +  e  1  J<1 

11  ,  .  a  ,  .  a  - 


1  +  4 


1  +  2 


EPU  (a) 


for  a  e  [0,1]  at.  ' 

.  fU  ,  .  ^~C01  ~Kl~hx2-X3  .  Xl+X3  x!+x2+x3 

min  fQ2VX,a)  5  e  +  c_2  e  +  e 

xeE3  1+1  (a- 2)  02 


subject  to  ("(x.a)  -  -4^-  e'2*1'2”2  +  - LffiS -  <  j. 


1  +  j  (a-2) 


1  +  f  (a-2) 


EPU2(a) 


for  a  e  (2/5,  1]  (for  a  €  [0,  2/5]  we  define  f^2(x»a)  =  +°°  for 
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k=0,l  ,  all  x  £  E  ).  The  underestimating  problem  is  defined  for 
a  £  [0,1]  as 

I  c  -V's‘2-*3  .  Xl+X3  .  *1+VX3 

min  f0(x,o) - 5 -  e  +  c02  e  +  c03  6 

xeE3  U  1+|  (a-1) 


i  T  cn/  ~2x  -2x0  4  %x_-x_ 

subject  to  f  1  (x,a)  =  - = -  e  + - ; -  e  <1 

1  +  ^-  (a-1)  1  +  ^  (a-1) 


EPL(a) 

The  above  three  problems  have  convex  optimal  value  functions 

f*U(ot)  ,  f*U(a)  ,  and  f*L(a)  ,  respectively,  for  a  c  [0,1]  . 

Recalling  from  Section  2  that  for  a  =  0  ,  cQ  =  (40,  20,  20,  1/3, 

4/3)  ,  5q  -  6*<c0)  =  (1/5,  2/5,  2/5,  3/10,  3/5)  ,  xQ  =  x*(cQ)  =  (0,  0, 

log2)  ,  and  f*(0)  ■  v*(0)  =  100.0  ,  application  of  inequality  (4)  gives 
us  convex  lower  bounds  on  f*(a)  ,  a  £  [0,1]  ,  as  follows: 


f*(a)  >  v(6Q,a)  = 


r  3/8  i 

601 

r  3/4 

o 
_ _ 

r  1/2  i 

605 

L1  + 1  <«-»J 

U  +  i  «*-i>_ 

L.1  +  \  (“-!)_ 

100.0 

.2  ;  :  i.3 


5  v  r  i  v  < 

1  +  ^  a  1  +  4  a  (1  +  a)- 


r  --  i 

611 

< 

r  1 

i14 

r  i  i 

«15 

Ll  +  |  (a-l)J 

Li +  \ 

i 

Li  +  |  (ot-i)J 

_ 95.24 

>1/3  r  3  31/6 


(i+fa)  [l+ja]  (1  +  ar 


The  linear  lower  bounds  on  f*(a)  can  be  easily  computed  using  the 
formulas 
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LQ(a)  -  f*(0)  +  Vav(60,0)a  =  100.0  -  86.67a 

and 

L^a)  -  f*(l)  +  Vav(61,l)(a-1)  -  51.96  -  21.65(a-l)  =  73.61  -  21.65a 

In  order  to  obtain  an  upper  bound  we  solve  the  dual  of  the  problem 
EPU^(a)  for  a  »  1  ,  obtaining  *  (320/13,  20,  20,  4/15,  8/9)  , 
r(l)  =  52/149  ,  and  6^  -  S*0^)  -  (45/149,  52/149,  52/149,  59/298, 

59/149)  .  From  this  we  calculate  f*U(l)  =  v(6^,c^)  =  68.34  and 
xj  -  x*U(cx)  -  (log(.9) ,  0,  log(1.325))  .  Since  f*U(0)  -  f*(0)  =  100.0  , 
the  linear  upper  bound  on  f*(a)  for  a  e  [0,1]  can  be  obtained  from 
the  formula 

U(a)  -  (l-a)f*u(0)  +  af*U(l)  *  (l-a)lOO.O  +  a68.34 

-  100.0  -  31.66a 

A  better  convex  upper  bound  can  be  obtained  In  the  form  fp1(xU(a),a)  , 
xU(a)  »  (l-a)xg  +  ax^  ,  a  e  [0,1]  .  The  above  bounds  are  depicted  in 
Figure  4. 
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Figure  4. — Bounds  on  f*(a). 
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